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(57) Abstract: The invention concerns a iterative method and a decoder for decoding space-time codes in a communication system 
with multiple transmission and reception antennae. The invention strikes a compromise between techniques based on interference 
cancellation algorithms such as BLAST, which show faulty performance concerning error rate based on signal-to-noise ratio and 
techniques based on maximum likelihood algorithms which are optimal in terms of performance, but highly complex in implantation 
such as the sphere decoder. Therefor the method consists in using a first matrix product between the received signal (Y) and a shaping 
matrix (B l ), and a second matrix product between a subtraction matrix (D 1 ) and the vector of the estimated symbols(S 1 ' 1 ) during the 
preceding iteration. The estimated symbols during the current iteration are generated by a subtracter (9) receiving the results (r l - z l ) 
of the two matrix products. The role of the matrix D 1 is to subtract from the current information symbol S l the interference caused 
by the other information symbols. 
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En ce qui concerne les codes a deux lettres et autres abrevia- 
tions, se referer aux "Notes explicates relatives aux codes et 
abreviations" figurant au debut de chaque numero ordinaire de 
la Gazette du PCT. 



(57) Abrege : L'invention concerne un proc&te iteratif et un decodeur lin£aire pour decoder des codes spatio-temporels dans un 
systeme de communication multi-antennes a remission et a la reception. L'invention est un compromis entre les techniques basees 
sur des algorithmes a annulation d' interferences tels BLAST, qui ont de mauvaises performances en terme de taux d'erreurs fonction 
du rapport signal a bruit, et des techniques basees sur des algorithmes a maximum de vrai semblance qui sont optimaux en terme 
de performances, mais d*une grande complexity d' implantation tel que le decodeur par spheres. Pour ce faire on utilise un premier 
produit matriciel entre le signal recu (Y) et une matrice de mise en forme (B 1 ) , et un second produit matriciel entre une matrice de 
soustraction (D 1 ) et le vecteur des symboles estim6s (S 1 " 1 ) lors de T iteration pr6ce<iente. Les symboles estim6s lors de rit^ration en 
cours sont g6ner6s par un soustracteur (9) recevant les resultats (r 1 , z l ) des deux produits matriciels. Le role de la matrice D 1 est de 
soustraire au symbole d' information S l courant T interference due aux autres symboles d' information. 
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"Precede pour decoder des codes espace-temps lineaires dans un 
systeme de transmission sans fil multi-antennes r et decodeur 
mettant en oeuvre un tel procede . " 



La presente invention se rapporte a un procede pour decoder 
des codes espace-temps lineaires dans un systeme de transmission 
sans fil multi-antennes. Elle concerne egalement un decodeur 
10 mettant en oeuvre ce procede . 

L' invention trouve une application particulierement 
interessante dans le domaine de la transmission ou de la diffusion 
radio de donnees numeriques, ou de donnees analogiques 
echantillonnees, notamment dans le cas de transmission avec les 
15 mobiles ou bien de fagon plus generale, dans le cas de reseau sans 
fil locaux ou non. D'une maniere plus precise, 1' invention peut 
notamment s'appliquer lorsqu' on desire faire des transmissions sans 
fils a grand debit. Une premiere categorie d' application concerne 
les systemes de communication cellulaire avec les mobiles tels que 
20 l'UMTS par exemple. Une seconde categorie d f application concerne 
les reseaux locaux sans fils. Une troisieme categorie d' application 
est celle des futurs reseaux ad hoc. 

Les systemes de communication sans fils actuels et a venir 
requierent une grande qualite de transmission pour des debits en 
25 constante augmentation. En effet f ces systemes ont pour object if la 
mise en place de services incluant voie, donnees et video de grande 
qualite . 

II existe de nombreuses contraintes liees a la propagation 
multi-chemins . Cette propagation multiple est due au faite qu f un 
30 signal radio atteint une antenne par plusieurs chemins via des 
reflexions sur 1' ionosphere et sur des objets terrestres tels que 
des montagnes et des batiments par exemple. Les effets de cette 
propagation multiple sont des interferences constructives et 
destructives ainsi que des decalages de phases du signal emis. 
35 Pour elaborer un systeme de communication fiable et robuste, 

on utilise des techniques de diversites spatiales, temporelles et 
f requentielles . La diversite spatiale consiste notamment en une 
disposition de plusieurs antennes a l f emission et/ou a la 
reception. Lorsqu' on y ajoute la diversite temporelle, on cree une 
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diversity spatio-temporelle n£cessitant un codage spatio-temporel k 
1' Emission ainsi qu'un decodage spatio-temporel a la reception. On 
sait que les codes espace-temps permettent d' augmenter les debits 
d f information transmis pour une qualite de service egale . 
5 Des codes spatio-temporels lineaires a grand rendement ont 

deja ete decodes par des algorithmes de decodage que 1' on peut 
classer en deux families : 

- les algorithmes a annulation d' interferences tel BLAST 
("Bell Laboratories Layered Space-Time" ) qui ont de 
10 mauvaises performances en terme de taux d' erreurs f onction 

du rapport signal a bruit. 

les algorithmes k maximum de vraisemblance qui sont 
optimaux en terme de performances, mais d'une grande 
complexite d' implantation tel que le d6codeur par spheres. 
15 On connait le document US 6 178 196 (SESHADRI et al) dans 

lequel est d§crit un systeme multi-antennes combinant le decodage 
a annulation d f interference avec le decodage a maximum de 
vraisemblance. Le systeme decrit a pour but d' isoler les signaux 
provenant d'une station donn6e en supposant que les symboles des 
20 signaux emis par une autre station sont correctement determines 
lors d'une estimation preliminaire. 

Un decodeur optimal pour les codes spatio-temporels 
lineaire a grand rendement a et6 trouv6 par Damen, Chkeif et 
Belfiore [O. Damen, A. Chkeif and J.-C. Belfiore, "Lattice Code 
25 Decoder for Space-Time Codes," IEEE Communications Letters, May 
2000]. Ses performances depassaient ce que l'on pouvait avoir 
avec le decodeur "BLAST" [G. D. Golden, G. J. Foschini, R. A. 
Valenzuela, P. W. Wolniansky, "Detection Algorithm and Initial 
Laboratory Results using the V-BLAST Space-Time Communication 
30 Architecture," IEE Electronics Letters, Jan. 7 1999]. Cependant, 
il souffre de trois probl£mes etroitement li£s k son 
implementation : 

1. II est d'une complexity assez importante, ce qui rend son 
implementation & grand debit tres probl£matique, 
35 2. Sa complexite depend fortement du rayon de la sphere choisie. 
Le mot de code decode doit en effet se trouver dans cette sphere 
si on veut que le decodage n' 6choue pas, et de plus, il faut que 
dans cette m§me sphere, on ne trouve pas beaucoup d f autres mots 
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de code, car sinon, la complexity de decodage peut encore 
augmenter de fagon tres iiaportante. 

3. Finalement, la complexity de decodage depend tres fortement du 
rapport signal a bruit et le bon choix de la sphere aussi. Les 
5 algorithmes de choix de la sphere sont eux-memes assez complexes. 
On voit done qu' il f allait trouver un autre algorithme de 
decodage qui ne souffre pas de ce genre de problemes. 

La presente invention a pour but 1' elaboration d' un nouveau 
decodeur spatio-temporel lineaire realisant un compromis entre la 
10 technique a annulation d' interference et la technique a maximum de 
vraisemblance. L' invention vise des performances tres proches de la 
technique a maximum de vraisemblance avec une implantation 
simplifiee par rapport a celle des algorithmes a annulation 
d r interference. 

On atteint les objectifs precitSs avec un procede iteratif 
pour decoder un ensemble de N signaux echantillonnes dans un 
systSme de communication spatio-temporel a M antennes d' emission 
et N antennes de reception- Les N signaux sont traites par 
- intervalles de temps T correspondant A la longueur temporelle du 
20 code spatio-temporel lineaire associe aux signaux d* emission. 
Dans chaque intervalle de temps T, les signaux echantillonnes 
regus sont representees sous forme d'une matrice signal Y. 



15 



peut etre mis sous la forme Y « H X + W (1) 
25 dans lequel X est le mot de code espace-temps emis et est 
represents par une matrice rectangulaire a M lignes et T 
colonnes, H est la matrice du canal et est representee par une 
matrice rectangulaire £ N lignes et M colonnes, et enfin W est le 
bruit, represents par une matrice rectangulaire & N lignes et T 
30 colonnes. 

On voit done que les changements de lignes correspondent a 
des changements d f antennes alors que les changement de colonnes 
correspondent a des changements d' instants d r <§ chant illonnage . 
Dans ce modele, les coefficients de la matrice X dependent 
35 lineairement des symboles d r information & transmettre, e'est-a- 
dire que les coefficients de la matrice transmise, x i3 avec i 
allant de 1 a M et j allant de 1 a T s'ecrivent sous la forme 
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p=l g=l 

avec s pq etant les symboles d r information et gij Pq etant des 
coefficients qui dependent du code choisi. Ainsi, tout code 
espace temps lineaire, c' est-a-dire tel que les mots emis X ont 

5 leurs coefficients Xij qui peuvent etre mis sous la forme 
prec£dente, peut etre decode par le precede selon l r invention. 

On va maintenant decrire un modele vectorise sur lequel le 
procede selon 1' invention pourra etre applique. Au lieu de 
travailler avec des matrices, on peut vectoriser 1' expression du 

10 signal regu dans 1' equation (1) et ecrire 

Y = vec (Y) = H X + W (2) 

Avec X = vec (X) = G S , G est de dimension (MT,MT) 

15 



H est une matrice (NT / MT) ayant sur la diagonale T fois la 
matrice H, les autres coefficients etant nuls. 

L' Equation (2) est alors equivalente k l f equation (1) et en 

20 constitue la version vectorisee. C'est sur cette version que le 

decodeur de 1' invention s f appliquera. II suffit maintenant de 
def inir la matrice du canal etendu 



25 

que l'on utilisera dans la suite. On peut done reecrire la 
version vectorisee du signal regu, 

Y = C S + W 

30 

oti Y est un vecteur colonne a NT composantes, C est une matrice 
NTxMTf S est le vecteur colonne des symboles a MT composantes et 
W est le vecteur colonne bruit a NT composantes. 



H 0 



Et 



0 H 
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= H G 
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Le procede iteratif de decodage selon la presente invention 
permet d'obtenir une estimation des symboles des signaux ends. Ce 
procede s' applique en particulier pour N superieur ou egal a M. 
Selon 1' invention, il comprend les etapes suivantes : 
5 . Pre-traitement du vecteur Y pour maximiser le rapport 

signal sur bruit+interf erence de facon a obtenir un signal r 1 , 

soustraction du signal r l par un signal z l au moyen 

d'un soustracteur, le signal z< etant obtenu par un post- 
traitement de reconstruction de 1' interference entre symboles a 
10 partir des symboles estimes lors de 1' iteration precedente, 

detection du signal genere par le soustracteur de 
facon a obtenir, pour 1' iteration en cours, 1' estimation des 
symboles des signaux emis. 

L'etape de pre-traitement peut etre realisee en operant une 
15 multiplication matricielle entre le vecteur signal Y et une matrice 
B, la matrice B etant mise a jour a chaque iteration. 

L'etape de post-traitement peut egalement etre realisee en 
operant une multiplication matricielle entre le vecteur des 
symboles estimes lors de 1' iteration precedente et une matrice D, 
20 la matrice D etant mise a jour a chaque iteration. 

Selon 1' invention, pour chaque iteration, un coefficient de 
correlation normalise p est calcule, la mise a jour d'une matrice 
etant obtenu en determinant de nouveaux coefficients de la matrice 
en fonction du coefficient de correlation obtenu pour 1' iteration 

25 precedente . 

Suivant un mode de realisation pr6fer6 de 1' invention, les 
N signaux sont traites par intervalles de temps T correspondant a 
la longueur temporelle du code spatio-temporel lineaire associe 
aux signaux d' emission, et l'etape de pre-traitement fait 
30 intervenir la matrice B pour maximiser le rapport signal sur 
bruit+interference dont la fonction de transfert est 
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avec V 1 — 



-i-i 



El 

E s 



; j£ = diag(c H V-C) ; 



I : indice d f iteration; p : coefficient de correlation normalise 
entre les symboles reels et les symboles estim6s; N 0 : variance 
du bruit; Es : energie moyenne d'un symbole; C : matrice du canal 
etendu . 

De meme, l'etape de post-traitement peut faire intervenir 
une matrice D pour la reconstruction de 1 1 interference entre 
symboles dont la fonction de transfert est : 



C-p^-Diag 



No 



]<i<MT 



20 



25 



30 



Avant la premiere iteration, on a aucune information sur 
les symboles. La matrice B a pour role de maximiser le rapport 
signal/ (bruit + interference). Le vecteur z est vide. Une 
representation du vecteur des symboles S peut etre telle que le k 
ieme symbole est : 

S k (regu) - S k (emis) + Y*CliSk &mis ) + bruU 

A la deuxieme iteration, la matrice B maximise tou jours le 
rapport SINR. La matrice D va mimer 1 1 interference entre symbole, 
c'est a dire *2uCliS fe mts ) lors de 1' iteration precedente au 
niveau du signal r sortant de la matrice B. Le soustracteur 
permet de retrancher cette interference. 

A la derniere iteration, on suppose que les symboles sont 
correctement estimes, c'est k dire que D permet de reconstruire 
toute 1' interference, de sorte qu'on estime que : 
S k (regu) = S k (emis) + bruit 

L' invention est un decodeur pouvant s 1 adapter aux codes 
espace-temps lineaires, quels qu'ils soient, c'est-£-dire qu' elle 
permet de decoder tout code tel que les sequences 
s f ecrivent comme une combinaison lineaire des 
d' information. 



emises 
symboles 
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On peut remarquer que les matrices B et D dependent du 
coefficient de correlation qui est different pour chaque 
iteration. Selon une caracteristique avantageuse de 1* invention, 
pour determiner le coefficient de correlation p* , a chaque 
iteration: 

- on calcule le rapport signal sur interference SINR a 
1' entree du detecteur a seuil en utilisant la formule suivante : 



SINR e = 



1 



-1 



avec 



10 



et 



1-pl-l 

1 • exponentiel integrale -^1 { s ) = £ ~f < ^ 

No 



NE, 



- on calcule la probabilite d'erreur symbole Pr, par 
exemple a 1' entree du detecteur a seuil, a partir du rapport 
15 signal sur interference SINR 4 ; et 

on calcule alors le coefficient de correlation p* a partir 
de la probabilite d'erreur symbole Pr. Pour ce faire, on peut 
utiliser une formule conventionnelle donnant le coefficient de 
correlation en fonction de la probabilite d'erreur symbole, cette 
20 formule dependant de la modulation utilisee a 1' Emission. 

On approxime ainsi de fagon la plus precise possible le 
coefficient de correlation p e . 

A la suite d' experimentations, quatre iterations ont suffit 
pour obtenir de tres bons r6sultats. Mais on peut aussi definir 
25 une valeur minimum (valeur seuil) du coefficient de correlation 
pour laquelle on interromps les iterations. Ce coefficient est 
essentiellement fonction de H, qui peut etre estime de facon 
conventionnelle, et de la variance N 0 . Par consequent, toutes les 
valeurs du coefficient de correlation p, ainsi que les valeurs 
30 des matrices B et D (pour 1' ensemble des iterations) peuvent etre 
calcuiees avant la premiere iteration. On peut done les stocker 
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avant la premiere iteration dans une memoire de facon a les 
ressortir ensuite pour chaque iteration. 
De preference, on pose p° =0. 

Par ailleurs, pour calculer la probability d'erreur symbole 
5 Pr, on peut supposer que le bruit total est gaussien et on peut 
utiliser la formule correspondant a la constellation d'une 

modulation lineaire, par exemple fi^^ P our le s y st ^e 

d'encodage BPSK ("Binary phase Shift Keying") ou bien utiliser 
des tables indiquant la probabilite d'erreur en fonction du 
10 rapport signal sur bruit. En effet, selon la complexity de la 
probabilite d'erreur par symbole il peut etre interessant de 
tabuler directement la formule. 

Avantageusement, pour calculer le coefficient de 
correlation p* a partir de la probabilite d'erreur symbole Pr, on 
15 peut supposer que lorsqu'il y a erreur, le detecteur a seuils 
detecte l'un des plus proches voisins du symbole emis. 

A titre d' exemple, a 1* iteration finale, on peut introduire 
le signal sortant du soustracteur dans un decodeur a entrees 
souples . 

20 Suivant un mode de mise en oeuvre pref^re de 1' invention, 

les symboles d' information peuvent etre des elements d'une 
constellation . issue d'une modulation d' amplitude en quadrature 
MAQ, ou QAM en langue anglaise. 

Suivant un autre aspect de 1' invention, il est propose un 
25 decodeur spatio-temporel pour decoder un vecteur signal Y obtenu 
a partir de N signaux echantillonnes dans un systeme de 
communication spatio-temporel a M antennes d' emission et N 
antennes de reception, avec N superieur ou egal a M, en vue 
d'obtenir une estimation des symboles des signaux emis. Selon 
30 1' invention, ce d6codeur comprend : 

un module de pre-traitement du vecteur Y pour 
maximiser le rapport signal sur bruit+interference de facon a 

obtenir un signal , 

un soustracteur pour soustraire un signal z l du signal 



35 r l 
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un module de post-traitement pour la reconstruction de 
1' interference entre syniboles a partir des symboles estimes lors 
de 1 T iteration precedente de fagon a generer le signal z l , 

un detecteur a seuil pour la detection du signal 
5 genere par le soustracteur de fagon a obtenir, pour 1 f iteration 
en cours, 1 1 estimation des symboles des signaux emis. 

Ces modules de pr6-traitement et de post-traitement peuvent 
etre des matrices, B et D, selon les formules indiquees 
precedemment . 

10 D'autres avantages et caracteristiques de 1' invention 

apparaitront £ l f examen de la description detaillee d' un mode de 
mise en ceuvre nullement limitatif, et des dessins annexes sur 
lesquels : 

- La figure 1 est un schema illustrant quelques Elements d'une 
15 chaine de transmission au sein d r un emetteur et d' un recepteur, 

le decodeur spatio-temporel selon 1' invention etant integre dans 
le recepteur; 

- La figure 2 est un schema general illustrant l f architecture du 
decodeur spatio-temporel selon 1' invention; et 

20 - La figure 3 est un schema bloc general du decodeur selon 
1» invention. 

Sur la figure 1 est represents un emetteur 1 dote d f une 
plurality d'antennes 7. L' 6metteur 1 comprend notamment en amont 
des antennes, un codeur correcteur d r erreurs 3 suivi d' un codeur 
25 spatio-temporel lineaire 4. Selon une caracteristique avantageuse 
de l f invention, le decodeur selon 1' invention est apte a 
s'appliquer & tout code espace-temps lineaire, c r est a dire des 
codes tels que les sequences 6mises s'ecrivent comme une 
combinaison lineaire des symboles d f information. 
30 Les signaux emis par les antennes 7 sont captes par une 

pluralite d'antennes 8 au sein d'un recepteur 2. Les signaux regus 
subissent un traitement au sein d'un decodeur spatio-temporel 5 
selon 1' invention de fagon a estimer les symboles d' information des 
signaux 6mis. Le decodeur spatio-temporel 5 prSsente une sortie 
35 souple a laquelle est connects un decodeur a entree souple pour 
decoder des codes correcteurs d' erreurs tels que des codes 
convolutifs, turbo-codes, code de Reed-Solomon,... ; le decodage 




WO 03/107582 PCT/FR03/01783 



- 10- 

pouvant se faire par un algorithme de Viterbi, un algorithme MAP 
(maximum a posteriori) ou LOG-MAP itere,.... 

Sur la figure 2 est illustre un schema general de 
1' architecture du decodeur spatio-temporel selon 1' invention. Ce 
5 decodeur met en oeuvre un procede iteratif permettant de determiner 
les symboles d' information S a la suite d'une decision ferme 
generee par un detecteur a seuil 10. Toutefois, le decodeur spatio- 
temporel selon 1' invention genere egalement des symboles 
d' information a la suite d'une decision souple apte a etre injecte 
10 dans le decodeur 6 a entrees souples, le signal injecte etant le 
signal obtenu lors de la derniere iteration. L' architecture du 
decodeur selon 1' invention fait principalement intervenir deux 
modules B et D injectant leurs signaux de sortie dans un 
soustracteur 9. Le soustracteur 9 genere un vecteur de symbole 
15 d' information a la suite d'une decision souple, ces symboles etant 
ensuite detectes par le detecteur a seuil 10 afin d' obtenir des 
symboles estimes par decision ferme. 

Les deux modules B et D representent des produits matriciels 
de leurs signaux d' entree par les matrices B« et D 4 , l'indice I 
20 traduisant 1' iteration en cours . Le module B recoit en entree le 
signal Y provenant des antennes 8. Le produit de la matrice BL par 
Y est un signal r l auquel on soustrait le signal z' issu du 
produit matriciel de la matrice D< par le signal ST 1 . Les matrices 
B e et sont telles que : 
25 en reprenant les notations de la figure 2, il vient alors : 

r*=B e Y 

En remarquant que dans la detection du symbole Sl du vecteur 
30 s, la i i6me composante de z' ne doit pas etre soustraite de r l , 
on impose la contrainte suivante sur D* : 



Z^=0, \/\<i<n 



\ 
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c'est a dire une diagonale nulle pour ne pas ret rancher le 
symbole utile. 

Le detecteur genere la decision dure S e 4 partir de S C e 
qui se traduit a 1' iteration i par : 



s l = B e Y-D* -s<- 1 =B e {H-s +w)-D e -s e - 1 



La premiere etape a mettre en oeuvre dans le decodeur itere 
consiste a determiner les matrices B< et D< telle que l'erreur 
10 quadratique moyenne a 1' entree du detecteur a seuil soit la plus 
petite possible. Elle est definie a 1' iteration 1 par la 
grandeur : 

MSE'iB'tD^^s' -s\ 2 

15 Minimiser 1' expression suivante implique que B< et D« verifient 

dMSEtjB',!)*) _ dMSE e (B',D e ) _ 
m dD* 

La resolution de ces equations permet d'obtenir, a 
20 1' iteration I en fonction de 1' iteration £-1 des matrices B« et I* 
suivantes : 



1 



N 0 



\<d<MT 



.C»V l 



25 tf^B'.Cp^-Diag 



_L 



&s J 



\<d<MT 



Avec : 
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et 



-1-1 



h£^±c-C H +Id N 



No. 
E s 



A* = diag(c H V<-C) 

ou e s est l'energie moyenne de la constellation MAQ et N 0 est la 
variance du bruit. 

La forme de la matrice & est intuitivement satisfaisante . 



En effet, si S* *=s tel que 



10 



15 



Ps 1= =l, alors la sortie de D* 
reproduit la composante intersymbole pour chaque symbole s if 

£-1 -I 

D'une maniere plus generale, P ~ A indique la confiance 
que l'on a dans la qualite de 1' estime s tA . Si n'est pas 

t—l i 

fiable alors P = A sera bas et par consequent une ponderation 
plus faible sera appliquee a 1' estimateur de 1' interference inter 
symbole soustraite de r l . A contrario, si 5 <_l est un excellent 
estime de s, alors p^-^let presque toute 1' interference inter 
symbole est soustraite de r' . A remarquer egalement qu' au 

premier passage e=l, Ps" 1 = PS = 0 et comme il n' y a pas d' estime 
s l ~ x =5° encore disponible, B« est le filtre adapte. 

Dans ces equations, p* designe la valeur, a 1' iteration I, de 
la correlation entre les symboles detectes a 1' iteration I, si et. 
les symboles reellement emis. Cette correlation est done 

;_ Ejs k sj) 
H E s 

ou E s est l'energie moyenne d f un symbole 
25 On peut remarquer que les matrices B* et D* ont des roles 

tres diff^rents : 



20 




WO 03/107582 -13- PCT/FR03/01783 

La matrice B e est une matrice de mise en forme du 
signal. A l r iteration zero, on constate que B e est le decodeur 
lineaire qui minimise l'erreur quadratique moyenne. Lorsque la 
correlation devient grande (tend vers 1), le role de la matrice 
5 B e devient marginal. 

Le role de la matrice D« est de soustraire au symbole 
d' information courant 1' interference due aux autres symboles 
d f information. Son role est marginal & 1' iteration z6ro, mais a 
mesure que la confiance sur les symboles detectes S e crolt, son 
10 role devient determinant. 

On va maintenant d£crire un moyen d f estimer le coefficient 
de correlation. 

h 

rapport signal sur interferences a 1' iteration I, SINR*. Ce 
15 rapport signal sur interferences ("Signal to Interference Noise 
Ratio" en langue anglaise) est d'abord calcule en fonction de la 
correlation a 1' etape precedente £-1. 

1. Calcul de SINR« 
On montre que : 



20 



Le calcul de la corr61ation p'=— ^ — reclame le calcul du 



r \ 

SINR 1 = 



V 



avec 



N 



1 l-{p t l Y 1 = N-E S _ 

Ei (s)= — —dt f i' exponentielle integrals 



25 et **' ' 1 



2. Calcul de 



Le calcul de f/ se fait en plusieurs etapes 

30 
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a- Poser p° =0 et £=1 

b- Calculer le rapport signal sur interferences (+ bruit) 
l 1 entree du detecteur a seuil compte tenu de la formule : 

5 



c- Calculer la probability d' erreur symbole Pr a 1' entree du 
10 detecteur a seuil a partir du SINR* en supposant que le bruit 
total est gaussien et en utilisant la formule correspondant a la 
constellation, par exemple dans le cas d'une modulation de type 
m-PSK, on peut utiliser la formule suivante : 



15 Ce calcul resulte d'une approximation a fort rapport signal a 
bruit pour les probabilites d' erreur symbole des modulation m-PSK 
associees a un detecteur a seuil symbole par symbole en presence 
d'un canal additif de type bruit blanc gaussien. Cette 
approximation est liee a la definition du modele equivalent. 

20 

d- Calculer 1' expression de f/ a la sortie du soustracteur en 
supposant que lorsqu' il y a erreur, alors le detecteur a seuils 
detecte l'un des plus proches voisins du symbole emis . Dans le 
cas d' une m-PSK on obtient : 



e- Incrementer £ = £+1) et retourner a l'etape b. 

Comme on le voit sur la figure 2, 1' architecture du decodeur 
selon 1' invention fait apparaitre des decisions souples a la sortie 
du soustracteur 9. Ces decisions sont exploitees a 1' iteration 
finale. Elles peuvent etre injectees dans un decodeur a entrees 
souples 6. 
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En notant L 1' iteration finale, on peut approximer $ par 
1' equation suivante : 

"g L = K S + bruit 

5 

Avec K constante > 0 connue, S le vecteur colonne des symboles a 

MT composantes . 

Si le decodeur a entrees souples ou ponderees est un 
algorithme de Viterbi, il n' aura qu' a minimiser la distance 
10 euclidienne entre J* et KS sur tous les mots du code correcteurs 
d' erreurs . 

Sur la figure 3 est represents un schema bloc general du 
decodeur selon 1' invention. Le module de decodage 11 comprend les 
elements B, D, 9 et 10 de la figure 2. On retrouve en effet le 
15 vecteur Y en entree, et une sortie souple g (sortie du 



20 



soustracteur 9) ainsi que la sortie s e (sortie du detecteur a 
seuil 10) . On distingue un module 12 de calcule des coefficients 
des matrices B et D. Ce module 12 permet de calculer, a chaque 
iteration ou de facon global avant le debut des iterations, le 
coefficient de correlation et les coefficients b£ et de des 
matrices B et D. Ce module 12 recoit en entree la matrice H 
(estimation du canal de transmission) et la variance N 0 . H peut 
generer en sortie le coefficient de correlation que 1 ' on peut 
utiliser pour interrompre les iterations ou pour toute autre 

25 utilisation. 

Bien sur, 1' invention n'est pas limitee aux exemples qui 
viennent d'etre decrits et de nombreux amenagements peuvent etre 
apportes a ces exemples sans sortir du cadre de 1" invention. 
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REVEND I CAT I ONS 



10 



15 



20 
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1. Procede itgratif pour decoder un vecteur signal Y obtenu §. 
partir de N signaux 6chantillonnes dans un systdme de 
communication spatio-temporel a M antennes d' emission et N 
antennes de reception, avec N supSrieur ou egal & M, en vue 
d'obtenir une estimation des symboles des signaux emis; 
caract£rise en ce que chaque iteration comprend les Stapes 
suivantes : 

Pr£- trait ement du vecteur Y pour maximiser le rapport 
signal sur bruit+interf Srence de fagon §. obtenir un signal r 1 , 

soustraction du signal r e par un signal z l au moyen 

d ! un soustracteur, le signal z l Stant obtenu par un post- 
traitement de reconstruction de 1 1 interference entre symboles k 
partir des symboles estimes lors de 1" iteration precedente, 

detection du signal g6nSr6 par le soustracteur de 
fagon a obtenir, pour 1 ! iteration en cours, l 1 estimation des 
symboles des signaux Smis; 

et en ce que, les N signaux etant trait^s par intervalles de 
temps T correspondant a la longueur temporelle du code spatio- 
temporel lineaire associe aux signaux d 1 emission/ l'etape de pr6- 
traitement fait intervenir une matrice B pour maximiser le 
rapport signal sur bruit+interf erence dont la fonction de 



transfert est : jg £ =Diag 



1 



As J 



\<d<MT 



,C H V £ 



avec V*= 



u:c H +id N 



; j[ =diag(C"-V'-C) ; 



I : indice d' iteration; p : coefficient de correlation normalise 
entre les symboles r£els et les symboles estimes; N 0 : variance 
du bruit; Es : energie moyenne d'un symbole; C : matrice du canal 
Stendu; 
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et en ce que 1 1 6tape de post-traitement fait intervenir une 
matrice D pour la reconstruction de 1 ' interference entre symboles 
dont la fonction de transf ert est : 

/ \ 



10 



D^B'.C-p^-Diag 



1 



&s J 



\<i<MT 



2. Precede selon la revendication 1, caracterise en ce que l ! etape 
de pre- trait ement est realis^e en operant une multiplication 
matricielle entre le vecteur signal Y et la matrice B, la matrice B 
etant mise k jour & chaque iteration. 



3. Procede selon la revendication 1 ou 2 , caracterise en ce que 
l'etape de post-traitement est r^alisee en operant une 
multiplication matricielle entre le vecteur des symboles estimes 
lors de 1' iteration pr^cSdente et la matrice D, la matrice D etant 

15 mise a jour a chaque iteration. 

4. Procede selon la revendication 2 ou 3, caracterise en ce que 
pour chaque iteration, le coefficient de correlation normalise p 
est calcuie, la mise a jour d ! une matrice etant obtenu en 

20 determinant de nouveaux coefficients de la matrice en fonction du 
coefficient de correlation obtenu pour 1 • iteration precedente. 

5. Procede selon l'une quelconque des revendications precedentes, 
caracterise en ce que pour determiner le coefficient de 

25 correlation pour chaque iteration : 

on calcule le rapport signal sur interference SINR en 

\ 



utilisant la formule suivante : $INR I gig? Jj^ Jjj*^ 



± 



et en 



def inissant 1 1 exponentiel integrale 
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on calcule la probabilite d'erreur symbole Pr k partir 
du rapport signal sur interference SINR* ; et 

on calcule alors le coefficient de correlation §. 
5 partir de la probabilite d'erreur symbole Pr. 

6. Procede selon la revendication 5, caracterise en ce qu f on pose 
p° =0. 

10 7. Procede selon la revendication 5 ou 6, caracterise en ce que 
pour calculer la probability d'erreur symbole Pr on suppose que 
le bruit total est gaussien. 

8. Precede selon la revendication 7, caracterise en ce qu'on 
15 utilise la formule correspondant a la constellation issue d'une 

modulation lineaire a l 1 emission. 

9. Procede selon l'une quelconque des revendi cat ions 5 a 8, 
caracterise en ce que, pour calculer le coefficient de 

20 correlation p* h partir de la probabilite d'erreur symbole Pr, on 
suppose que lorsqu'il y a erreur, le detecteur & seuils detecte 
l'un des plus proches voisins du symbole emis. 

10. Procede selon l'une quelconque des revendi cat ions 
25 precedentes, caracterise en ce qu'& I 1 iteration finale, on 

introduit le signal sortant du soustracteur dans un decodeur a 
entries souples . 

11. Procede selon l'une quelconque des revendications 
30 pr6cedent es # caracterise en ce que les symboles d' information 

sont des elements d'une constellation issue d'une modulation 
d' amplitude en quadrature. 

12. Decodeur spatio-temporel mettant en oeuvre un procede selon 
35 l'une quelconque des revendications precedentes, pour decoder un 

vecteur signal Y obtenu & partir de N signaux 6chantillonnes dans 
un systeme de communication spatio-temporel a M antennes 
d' emission et N antennes de reception, avec N superieur ou 6gal & 
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M, en vue d'obtenir une estimation des symboles des signaux emis; 
caracterise en ce qu'il comprend : 

un module de pr£-traitement du vecteur Y pour 
maximiser le rapport signal sur bruit+interf erence de fa<?on a 

obtenir un signal r i , 

un soustracteur pour soustraire un signal z l du signal 

r l , 

un module de post-traitement pour la reconstruction de 
1 1 interference entre symboles a partir des symboles estim^s lors 
de I 1 iteration prec£dente de fagon k gen€rer le signal z l , 

un detecteur seuil pour la detection du signal 
gener£ par le soustracteur de fagon obtenir, pour 1 1 iteration 
en cours, 1' estimation des symboles des signaux emis; 
et en ce que, les N signaux etant traites par intervalles de 
temps T correspondant la longueur temporelle du code spatio- 
temporel lineaire associe aux signaux demission, le module de 
pre-traitement consiste en une matrice B pour maximiser le 
rapport signal sur bruit+interf erence dont la fonction de 



transf ert est : £ =Diag 



A 2 -i4' + 



avec V l ~ 



El 
E s 



No. 
E s j 



]<i<MT 



.C»V e 



-l-i 



N 



; X =diag(C"-V e 'C) ; 



I indice d 1 iteration; p : coefficient de correlation normalise 
entre les symboles reels et les symboles estimes; N 0 : variance 
du bruit; Es : energie moyenne d'un symbole; C : matrice du canal 
etendu; 

et en ce que le module de post-traitement consiste en une matrice 
D pour la reconstruction de 1 ' interference entre symboles, dont 
la fonction de transfert est : 



\<i<MT 
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13. Decodeur selon la revendication 12, caracterise en ce qu'il 
comprend tin decodeur a entries souples recevant le signal issu du 
soustracteur lors de l 1 iteration finale. 
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Categorte • Identification des documents cites, avec, le cas 6cheant, Flndlcation des passages pertinents 



no. des revendlcatlons visees 



WO 01 19013 A (HOME WIRELESS NETWORKS INC) 

15 mars 2001 (2001-03-15) 

page 2, Hgne 29 - Hgne 32 

page 3, Hgne 6 - Hgne 22 

page 7, colonne de gauche, ligne 30 - 

ligne 37 

page 9, colonne de gauche, ligne 4 
-colonne de drolte, ligne 12 
page 12, colonne de gauche, ligne 32 - 
ligne 52; figure IB 

EP 1 069 722 A (HUSHES ELECTRONICS CORP) 
17 janvier 2001 (2001-01-17). 
page 11, alinea 56 - alinea 57 
page 14, alinea 60; figure 5 
page 15, ligne 66 
page 16, Hgne 70 
figure 5 

-/- 
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| y| voir la suite du cadre C pour la fin de la Bste des documents 



Les documents de families de brevets sont lndlques en annexe 



° Categories speclales de documents cites: 

■A' document defintesant Fetal general de la technique, non 

consldere comme partlculierement pertinent 
•E" document anterieur, mate publle a la date de dep6t International 

ou apres cette date 
'L' document pouvant Jeter un doute sur une revendlcatlon de 

priorite ou cite pour determiner la dale de publication tfune 

autre citation ou pour une raison speclale (telle qu'indiquee) 
■O" document se rSferant a une divulgation orale, a un usage, a 

une exposition ou tous autres moyens 
*P* document publle avant la date de depot Internationa], mals 

posterieurement a la date de prtortte revendiquee 



■T document ulterleur publle apres la date de d6p0t International ou la 
date de prtortte et n'appartenenant pas a retat de la 
technique pertinent, mals cite pour comprendre le princlpe 
ou la theorte constfiuant la base de flnventfon 

•X' document partlculierement pertinent; Plnven tlon revendiquee ne petit 
etre consld6ree comme nouvelle ou comme Impltauant une activlte 
Inventive par rapport au document considers Isolement 

•V document partlculierement pertinent; Plnven tlon revendtauee 
ne peut etre conslderee comme Impllquant une activlte Inventive 
lorsquele document est assocle a un ou plusleurs autres 
documents de meme nature, cette comblnaison etant evldente 
pour une personne du metier 

document qui fait parte de la meme famllle de brevets 



Date a laquette la recherche Internationale a ete effectlvement achevee 



6 janvier 2004 



I Nom et adresse postals de Fadmlnlstratlon chargee de la recherche Internationale 
Office Europeen des Brevets, P.B. 581 8 Patentlaan 2 
NL-2280 HV Rljswljk 
Tel. (+31-70) 340-2040, Tx, 31 651 epo nl, 
Fax (+31-70) 340-3016 



Date d'expedltlon du present rapport de recherche Internationale 

14/01/2004 



Fonctlonnalre autorise 



Papantonlou, A 



FoCTUitalre PCT/1SA/21 0 (deuxifema teulDe) QuDlet 1992) 



HAPPORT DE REQMERCHE INTERNATIONALE 



C(sulte) DOCUMENTS CONSIDERES COMME PERTINENTS 



Categoric °j Identification des documents cites, avec.le cas echeant, rindlcatlondes passages pertinents 



Demand: 

PCT/FR 



fr m 



ale No 
1783 



no. des revendications vlsees 



BAUCH G ET AL: "Iterative equalization 
and decoding with channel shortening 
filters for space-time coded modulation" 
VEHICULAR TECHNOLOGY CONFERENCE FALL 2000, 
24 - 28 septembre 2000, pages 1575-1582, 
XP010524303 
Boston, MA, USA 

page 1576, colonne de gauche, llgne 1 - 
Hgne 9 

page 1576, colonne de droite, ligne 1 - 
Hgne 7 

page 1576, colonne de droite, ligne 20 
-page 1577, colonne de gauche, ligne 1 
page 1577, colonne de gauche, ligne 10 - 
ligne 24 

page 1580, colonne de droite, ligne 13 - 
ligne 17 

page 1580, colonne de droite, ligne 20 - 
ligne 29; figure 6 
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Foimulalre PCT/1SA/210 (suite de la deuxfeme teuilla) (iuIDet 1992) 



HAKKUKI Ufc HbUMtHUMt IN I tKIMAI lUNALfc 
Renselgnements relaUfs aux rnemo^^ftamllles de brevets 



PCT/FR 



fe n 



nale No 

1783 



Document brevet cite 
au rapport de recherche 



Date de 
publication 



Membre(s) de la 
famllle de brevet(s) 



Date de 
publication 



WO 0119013 



EP 1069722 



15-03-2001 



AU 

WO 



7358100 A 
0119013 Al 



10-04-2001 
15-03-2001 



17-01-2001 EP 



1069722 A2 



17-01-2001 



Formutaire PCT/ISA/210 (annexe families de brevets) (julDet 1992) 
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